Preconditioning CGNE iteration for inverse problems
نویسنده
چکیده
The conjugate gradient method applied to the normal equations (cgne) is known as one of the most efficient methods for the solution of (non-symmetric) linear equations. By stopping the iteration according to a discrepancy principle, cgne can be turned into a regularization method. We show that cgne can be accelerated by preconditioning in Hilbert scales, derive (optimal) convergence rates with respect to data noise for the preconditioned method, and illustrate the theoretical results by numerical tests.
منابع مشابه
A New Inexact Inverse Subspace Iteration for Generalized Eigenvalue Problems
In this paper, we represent an inexact inverse subspace iteration method for computing a few eigenpairs of the generalized eigenvalue problem Ax = Bx [Q. Ye and P. Zhang, Inexact inverse subspace iteration for generalized eigenvalue problems, Linear Algebra and its Application, 434 (2011) 1697-1715 ]. In particular, the linear convergence property of the inverse subspace iteration is preserved.
متن کاملPreconditioned Krylov subspace methods for the solution of least-squares problems
and Kk(BA,Br) = span{Br, (BA)Br, . . . , (BA)k−1Br}, (3) where B ∈ Rn×m is the mapping and preconditioning matrix, and apply Krylov subspace iteration methods on these subspaces. For overdetermined problems, applying the standard CG method to Kk(BA,Br) leads to the preconditioned CGLS [3] or CGNR [9] method while for underdetermined problems it leads to preconditioned CGNE [9] method. The GMRES...
متن کاملTR-2008006: Additive Preconditioning, Eigenspaces, and the Inverse Iteration
We incorporate our recent preconditioning techniques into the classical inverse power (Rayleigh quotient) iteration for computing matrix eigenvectors. Every loop of this iteration essentially amounts to solving an ill conditioned linear system of equations. Due to our modification we solve a well conditioned linear system instead. We prove that this modification preserves local quadratic conver...
متن کاملAdditive Preconditioning , Eigenspaces , and the Inverse Iteration ∗
We incorporate our recent preconditioning techniques into the classical inverse power (Rayleigh quotient) iteration for computing matrix eigenvectors. Every loop of this iteration essentially amounts to solving an ill conditioned linear system of equations. Due to our modification we solve a well conditioned linear system instead. We prove that this modification preserves local quadratic conver...
متن کاملGeneral constraint preconditioning iteration method for singular saddle-point problems
For the singular saddle-point problems with nonsymmetric positive definite (1, 1) block, we present a general constraint preconditioning (GCP) iteration method based on a singular constraint preconditioner. Using the properties of the Moore-Penrose inverse, the convergence properties of the GCP iteration method are studied. In particular, for each of the two different choices of the (1, 1) bloc...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Numerical Lin. Alg. with Applic.
دوره 14 شماره
صفحات -
تاریخ انتشار 2007